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Absence  of  magnetic  ordering  in  NiGa2S4 


1. 1.  Mazin 

Code  6390,  Naval  Reserach  Laboratory,  Washington,  DC  20375 
(Dated:  June  22,  2010) 

Triangular-layered  NiGa2S4,  contrary  to  intuitive  expectation,  does  not  form  a  noncollinear  an¬ 
tiferromagnetic  structure,  as  do  isoelectronic  NaCrCH  and  LiCrCH.  Instead,  the  local  magnetic 
moments  remain  disordered  down  to  the  lowest  measured  temperature.  To  get  more  insight  into 
this  phenomenon,  we  have  performed  first  principles  calculations  of  the  first,  second  end  third  neigh¬ 
bors  exchange  interactions,  and  found  that  the  second  neighbor  exchange  is  negligible,  while  the 
first  and  the  third  neighbor  exchanges  are  comparable  and  antiferromagnetic.  Both  are  rapidly 
suppressed  by  the  on-site  Hubbard  repulsion. 


NiGa2S4  occurs  in  the  layered  structure  where  the 
main  motif  is  a  triangular  layer  of  Ni2+  ions  surrounded 
by  edge-sharing  S  octahedra,  forming  a  trilayer  S-Ni-S 
with  the  rhombohedral  stacking  ABC,  and  the  trilayers 
are  separated  by  the  gallium  oxide  layers.  Ni2+  has  the 
electronic  configuration  f^e2,  and  one  expects  it  to  be 
insulating  and  magnetic  with  S  =  1,  and  the  magnetic 
moment  per  Ni  being  somewhat  less  than  2  ^b-  Indeed, 
this  is  exactly  what  happens  with  structurally  similar 
transition  metal  oxides  with  a  transition  metal  in  the 
ds  configuration,  such  as  NaCr02  or  LiCr02-  3d  metal 
ions  in  this  confuration  do  not  have  an  orbital  moment, 
therefore  one  expects  a  vanishing  single-site  anisotropy 
and  magnetic  interactions  reasonably  well  described  by 
the  Heisenberg  model.  In  the  nearest  neighbor  approx¬ 
imation  this  leads  to  noncollinear  ground  states,  with 
neighboring  spins  pointing  roughly  at  120°  to  each  other. 
Indeed,  this  is  what  has  been  observed  in  the  above- 
mentioned  chromates. 

NiGa2S4,  on  the  other  hand,  has  attracted  substan¬ 
tial  recent  interest  exactly  because  the  experiments  in¬ 
dicate  absense  of  any  long-range  magnetic  ordering^. 
Several  explanations  have  been  proposed,  such  as  full 
cancellation  of  the  nearest  neighbor  exchange  and  frus¬ 
trated  competition  between  the  2nd  and  the  3rd  neigh¬ 
bor  interactions^,  or  biquadratic  exchange^.  These,  how¬ 
ever,  impose  very  severe  quantitative  restriction  on  the 
exchange  parameters,  which  seem  quite  unrealistic. 

In  order  to  elucidate  magnetic  interactions  in  this  sys¬ 
tem  we  have  performed  ab  initio  density  functional  the¬ 
ory  (DFT)  calculations  of  the  electronic  structure  and 
magnetic  energies  of  NiGa2S4  and  found  that  indeed  con¬ 
dition  required  by  either  explanation  are  very  unlikely  to 
be  satisfied,  however,  the  magnetic  interactions  are  very 
rich  and  (as  conjectured  in  Refi)  long  range,  so  that  tak¬ 
ing  into  account  three  neighbor  shells  is  indispensable. 

For  the  calculations,  the  experimental  crystal 
structure^  was  used.  A  full  potential  linear  augmented 
plane  wave  code^  was  used  with  a  gradient  approxima¬ 
tion  for  exchange  and  correlation^.  The  calculations,  as 
expected,  render  an  insulating  band  structure  shown  in 
Fig.  1.  As  one  can  see,  Ni  d{eg)  bands  are  fully  po¬ 
larized,  a  small  gap  opens  (as  usual,  the  absolute  value 
of  the  gap  in  the  DFT  cannot  be  taken  very  seriously), 
the  magnetic  moment  per  Ni  is  2  /is,  and  in  the  calcu¬ 


FIG.  1:  LDA  band  structure  of  the  ferromagnetic  NiGa204. 
The  Ni  character  is  emphasized  by  the  size  of  the  circles. 


lations  this  moment,  not  unexpectedly,  resides  entirely 
in  the  MS2  layer.  Interestingly,  more  than  20%  of  the 
total  magnetic  moment  is  located  on  the  S  sites.  This 
creates  a  ferromagnetic  interaction  between  the  nearest 
neighbor  nickels,  which  in  the  theory  of  strongly  corre¬ 
lated  magnetic  systems  is  known  as  “ferromagneting  90° 
exchange” .  Note  that  in  the  DFT  the  Hund  rule  energy 
is  approximated  as  f  dr/(r)m2(r)/4  ss  ^  hMf  / 4,  where 
m( r)  is  the  total  spin  density,  I  and  Mj  are  the  Stoner 
factor  and  the  total  magnetization  of  the  atom  i ,  and  this 
energy  therefore  does  not  depend  on  the  Ni-S-Ni  bond  an¬ 
gle.  Let  us  estimate  this  interaction^.  Two  Ni  neighbors, 
when  their  spins  are  parallel,  induce  a  magnetic  moment 
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of  ~  0.2  /zb  on  each  of  the  two  bridging  sulfurs,  gaining 
an  additional  magnetic  energy  of  2/^0. 22/4  «  20  meV 
(The  Stoner  factor  of  the  sulfur  ion  can  be  estimated  as 
described  in  Ref0  and  is  about  1  eV).  Defining  J  as  half 
of  the  energy  for  flipping  a  bond  (Enn>  =  Jnn>Sn  •  Sn/), 
we  find  a  ferromagnetic  contribution  to  J  of  the  order  of 
10  meV.  This  is  not  a  small  energy,  but  it  has  to  compete 
with  also  large  conventional  antiferromagnetic  superex¬ 
change.  The  latter  in  DFT  is  2t2/Im,  where  t  is  the 
effective  Ni-Ni  d  —  d  hopping  and  I  Mi  is  the  Hund  en¬ 
ergy  cost  of  exciting  an  electron  with  the  opposite  spin, 
In  ~  0.8  eV.  However,  taking  into  account  on-site  Mott- 
Habbard  correlation  mandates  substiting  the  Stoner  I 
in  this  expression  by  the  Hubbard  U,  which  is  at  least 
4  times  larger.  While  Ni  in  NiGa2S4  is  not  necessarily 
strongly  correlated,  there  is  no  doubt  that  the  energy  of 
adding  an  electron  is  substantially  underestimated  in  the 
DFT,  and  hence  the  AFM  superexchange  overestimated. 


FIG.  2:  (color  online)  Different  magnetic  patterns  used  in 
calculating  the  exchange  constants.  Squares  (circles)  indicate 
up  (down)  moments  within  in  the  supercell,  thin  dark  (thick 
light)  line  antiferro-  (ferro-)magnetic  bonds.  The  first  pattern 
corresponds  to  the  \/3  x  \/3  supercell  indicated  by  the  dash 
lines,  the  second,  third  and  fourth  patterns  to  a  2  x  1,  3  x  1, 
and  2x2  supercells,  respectively 

Let  us  now  investigate  numerically  the  magnetic  inter¬ 
actions  in  NiGa2S4.  First,  we  want  to  make  sure  that  our 
conjecture  about  the  absence  of  magnetic  anisotropy  is 
indeed  correct.  This  can  be  adressed  by  running  fully 
relativistic  calculations  imposing  different  magnetic  field 
directions  and  comparing  energies.  The  result  is  that 
both  energies  differ  by  at  most  0.03  meV.  That  is  to 
say,  the  single  site  anisotropy  is  not  an  important  fac¬ 
tor  in  NiGa2S4.  Having  established  that,  we  have  com¬ 
puted  several  different  collinear  magnetic  structures,  as 
shown  in  Fig.  2.  If  mapped  onto  a  Heisenberg  model 
with  three  nearest  neighbor  interactions,  these  give  the 
exchange  constants  of  8.4,  0.3  and  4.1  meV  (defined  so 
that  the  total  energy  is  equal  to  sum  over  all  bands  of 
JijSiSji  |Sj|  =  1)  for  the  first,  second  and  third  neigh¬ 
bors,  respectively. 

Several  observations  are  in  place.  First,  in  LDA,  while 


FIG.  3:  (color  online)  Calculated  exchange  constants  for  the 
first  three  neighbor  shells  in  NiGa204,  in  meV,  as  a  function 
of  the  Hubbard  U,  assuming  the  intraatomic  J  =0.07  Ry.  The 
additional  entry  at  U  =  0  corresponds  to  J  =  0. 


the  second  neighbor  exchange  is  negligible,  the  third  one 
is  sizeable  and  comparable  with  the  nearest  neighbor 
exchange  (it  was  suggested  already  in  Ref.  Q]  that  the 
3rd  neighbor  exchange  may  be  anomalously  large  in  this 
compound).  This  can  be  traced  down  to  anomalously 
large  3rd  neighbor  hopping,  which  is  in  fact  generic  for 
triangular  layers  of  transition  metal  oxides.  Indeed,  if 
the  metal-oxygen  bonds  form  precisely  90°  angles,  the 
strongest  nearest  neighbor  hopping  channel,  eg  —  p  —  eg 
(or,  in  compounds  like  NaxCo02,  t2g  —  p  —  t2g)  is  for¬ 
bidden  by  symmetry,  however,  a  third  neighbor  path, 
eg  —  p  —  t2g  —  p  —  eg  is  fully  allowed  and  in  fact  has  the 
most  favorable  geometry  (Fig[4}8.  This  creates  a  possib- 
lity  for  a  sizeable  superexchange  of  the  order  of,  as  usu¬ 
ally,  t2jf/ A,  where  teff  is  the  effective  hopping  that  ap¬ 
pears  after  all  intermediate  states  are  integrated  out  and 
A  is  the  energy  required  to  flip  the  spin  of  a  metal  ion.  In 
LDA,  teff  is  of  the  order  of  t^dcr/ (Ed  -  Ep)/ (Eeg  -  Et2g), 
and  A  is  of  the  order  of  the  Stoner  (Hund)  parameter, 
<  1  eV.  In  the  Hubbard  model,  on  the  other  hand,  A  is 
set  by  the  scale  of  Hubbard  U  ~  4  —  6  eV.  As  usually, 
the  real  life  is  somewhat  in  between,  meaning  that  the 
exchange  constants  in  LDA  are  likely  overestimated. 

This  can  be  easily  demonstrated  using  the  LDA+U 
method  that  takes  into  account  the  Mott-Hubbard  cor¬ 
relations  in  some  crude  approximation  (Fig.  [3j) .  To  get 
an  idea  of  the  overall  scale  of  the  picture  we  have  esti¬ 
mated  the  value  of  U  using  the  quasiatomic  loop  in  a 
standard  Linear  Muffin  Tin  Orbital  package,  as  decribed 
in  Ref.  @,  and  obtained  U  ~  0.3 Ry.  This  simplified 
method  is  known  to  underestimate  U,  therefore  we  have 
carried  out  calculations  up  to  U  =  0.45  Ry.  At  that  max¬ 
imal  value  of  U  all  three  exchange  constants  practically 
vanish,  within  the  accuracy  of  the  calculation.  In  fact, 
the  nearest  neighbor  constant  at  U  =  0.45  Ry  becomes 
negative  (—0.04  meV),  but  for  all  practical  purpose  it 
may  be  considered  zero.  Interestingly,  at  this  value  of 
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FIG.  4:  (color  online)  Different  possible  superexchange  paths 
within  NiS2  trilayers.  Thin  dashed  triangles  indicate  the  tri¬ 
angular  Ni  and  (in  the  left  panel)  S  layers.  Thick  (red  online) 
lines  show  the  exchange  paths.  The  solid  parts  correspond  to 
short  bonds  (Ni-S,  2.422  Aand  the  shorter  S-S  bond,  3.212 
A),  the  dashed  ones  to  the  longer  S-S  bond,  3.626 
AA. 

U  the  sum  of  our  calculated  exchange  constants  over  all 
(six)  bonds  gives  1. 1*6=6. 6  meV=76  K,  to  be  compared 
with  the  Curie- Weiss  temperature  of  78±1  K—  (the  cal¬ 
culated  sign  is  antiferromagnetic,  in  agreement  with  the 
experiment) . 

What  prevents  the  system  from  ordering  remains  un¬ 
clear.  Nakatsuji  et  alX  conjectured  that  the  ratio  of 
the  nearest  and  the  third  neighbor  exchange  (no  second 
neighbors)  is  ~  —0.2,  while  our  ratio,  at  large  U ,  is  es¬ 
sentially  zero.  To  this  point  it  should  be  mentioned  that 
neither  the  accuracy  of  LDA+U  functional  is  sufficient 
to  make  firm  statements  with  a  precision  of  a  fraction 
of  a  meV,  not  the  three-shell  isotropic  Heisenberg  model 
is  accurate  to  that  extent.  If  the  calculated  numbers  at 
U  =  0.45  Ry  are  off  by  «  0.2  meV  this  would  be  enough 
to  bring  the  calculated  numbers  in  consistency  with  the 
Nakatsuji  et  aVs A  model. 

Finally,  one  should  keep  in  mind  that  the  actual  or¬ 
dering  temperature,  if  any,  is  suppressed  by  the  2D  char¬ 
acter  of  magnetic  interactions.  Indeed,  NiGa2C>4  has  an 
extra  GaC>2  layer  compared  to  typical  layer  oxides  of  the 
formula  ABO2  and  therefore  the  surexchange  interaction 
between  the  layers  is  reduced.  We  have  estimated  this 
interaction  by  comparing  the  energies  of  the  fully  fer¬ 
romagnetic  ordered  structrure  and  the  A-type  antiferro¬ 
magnetic  one  (FM  layers  stacked  antiferromagnetically) . 
The  former  is  higher  by  about  1  meV  (J  =0.5  meV)  in 
LDA  and  by  about  0.3  meV  (J  =  0.15  meV)  in  LDA+U 
{U  =  0.3  Ry).  Importantly,  this  superexchange  is  addi¬ 
tive  with  respect  to  all  possible  hopping  paths  between 
the  layers  that  include  not  only  hopping  from  a  Ni  to  the 
other  Ni  right  above,  but  also  to  a  large  number  of  neigh¬ 
boring  Ni  sites  in  the  next  planed.  The  real  exchange 
coupling  between  the  two  antiferromagnetic  planes  will 
be  additionally  reduced.  An  estimate  of  J±  ~  0.05  meV 
seems  reasonable.  Of  course,  although  a  computer  code 
calculates  the  numbers  with  arbitrary  precision,  the  ac¬ 
tual  physical  approximations  used  in  the  calculations  pre¬ 


clude  statements  about  energy  differences  of  the  order  of 
0.05  meV.  For  all  practical  purposes,  it  may  be  less  than 
0.01  meV,  which,  of  course  would  make  long-range  order¬ 
ing  at  the  experimentaly  probed  temperatures  (a  fraction 
of  a  Kelvin)  impossible. 

To  conclude,  we  have  calculated  the  three  near¬ 
est  neighbor  exchange  constants  in  NiGa204  by  map¬ 
ping  LDA+U  calculations  onto  the  isotropic  Heisenberg 
model.  We  found  that  (a)  an  anomalously  large  third 
neighbor  coupling  exists  in  the  system  that  can  be  traced 
down  to  superexchange  via  occupied  t2g  orbitals,  (b) 
upon  including  correlation  effects  in  terms  of  Hubbard 
U,  all  exchange  constants  decrease  rapidly,  according  to 
general  superexchange  intuition,  but  the  first  neighbor 
exchange  is  more  rapidly  suppressed  than  the  third  neigh¬ 
bor  one,  (c)  at  U  =  0.45  Ry  (6  eV)  the  nearest  neighbor 
exchange  is  entirely  suppressed  (and  we  could  not  ex¬ 
clude,  based  on  our  calculations,  that  it  does  not  become 
slightly  ferromagnetic),  while  the  third  one  is  exactly  the 
right  magnitude  to  explain  the  observed  Neel  temper¬ 
ature.  On  the  other  hand,  U  =  6  eV,  at  least  on  an 
intuitive  level,  seems  to  be  too  large  for  Ni  in  such  an 
enviroment.  At  more  realistic  C/’s,  such  as  3-4  eV,  both 
interactions  remains  firmly  antiferromagnetic.  It  should 
be  kept  in  mind  that  energy  differences  on  the  order  of  1 
meV  are  on  the  borderline  of  many  approximations  used 
in  our  analysis.  It  is  possible  that  several  weak  effects 
conspire  to  prevent  the  system  from  ordering.  First  of 
all,  in  particular  the  nearest  neighbor  interaction  results 
from  cancellation  of  two  considerably  stronger  superex¬ 
change  interactions  of  the  opposite  signs:  the  AFM  su¬ 
perexchange  and  the  FM  superexchange  due  to  the  Hund 
rule  coupling  on  S.  Albeit  we  see  no  obvious  reason  for 
the  DFT  to  underestimate  the  latter,  such  a  possibil¬ 
ity  cannot  be  excluded.  Second,  we  did  not  make  any 
attempt  to  evaluate  further  exchange  interactions  be¬ 
yond  the  third  shell.  While  there  is  no  special  mech¬ 
anism  making  them  sizeable  (as  opposed  to  the  third 
neighbor  interaction),  again  we  cannot  prove  that  nu¬ 
merically.  Third,  while  the  calculated  on-site  anisotropy 
is  very  small,  other  effects  beyond  the  isotropic  Heisen¬ 
berg  model,  such  as  biquadratic  exchange,  dipole-dipole, 
multispin  interactions,  eto^  while  small,  may  be  not  neg¬ 
ligible  on  the  backgound  of  the  strong  cancellation  of  the 
FM  and  the  AFM  superexchange,  and  may  possibly  cre¬ 
ate  additional  frustration  in  the  system.  Finally,  fourth, 
interplanar  coupling  between  the  antiferromagnetic  non- 
collinearly  correlated  planes  is  at  most  a  fraction  of  a 
Kelvin,  and  possibly  even  smaller.  This  additionally  su- 
presses  long  range  ordering. 

The  author  is  grateful  to  D.I.  Khomskii  for  useful  dis¬ 
cussion  and  to  C.  Broholm  for  bringing  up  this  problem. 
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corresponds  to  the  shorter  S-S  bond  and  t'pp  to  the  longer 
one.  Using  Andersen’s  and  Harrison’s  canonical  scalings, 
we  have  t'pp/tpp  «  (3.212/3.626)3  «  0.7,  and  t'ppcr/t'pp 17r  at  2. 
This  suggests  that  ts  at  t2  at  (2/3)fi.  This  is  obviously 
incorrect  and  indicates  that  (a)  the  hopping  parameters 
assumed  by  Takubo  et  al  do  not  agree  with  the  LDA  ones 
and  (b)  particularly  in  the  nearest  neighbor  hopping,  other, 
more  complex  paths,  such  as  Ni-S-Ga-S-Ni,  play  a  very  im¬ 
portant  role. 
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